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Critical behavior in short-time dynamics is investigated by a Monte Carlo study for the random- 
bond Potts ferromagnet with a trinary distribution of quenched disorders on two-dimensional tri- 
angular lattices. The universal dynamic scaling is verified and applied to estimate the critical 
exponents 6, z and j3/v for several realizations of the quenched disorder distribution. Our critical 
scaling analysis strongly indicates that the bond randomness influences the critical universality. 
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I. INTRODUCTION 

The study of the critical properties of physical sys- 
tems with bond randomness on phase transitions is a 
quite active field of current interest in equilibrium sta- 
tistical physics ||. One of the central importance 
here is to answer the question whether the critical ex- 
ponents (of a homogeneous pure magnet) change on the 
addition of quenched impurities and, if so, how do they 
change? Replying to this question, it was first stated al- 
ready two decades ago that if the specific heat critical 
exponent a of a pure system is positive, then a quenched 
disorder is a relevant perturbation at the second-order 
critical point and it causes changes in critical exponents. 
This statement is known as the Harris criterion 0. Fur- 
thermore, following the earlier work of Imry and Wor- 
tis Q] , who argued that quenched disorder could smooth 
first- order phase transitions and thus produce the second- 
order phase transitions, the introduction of randomness 
to pure systems undergoing a first-order transition has 
been comprehensively considered j^, [|. The theory was 
initially numerical checked with the Monte Carlo (MC) 
method by Chen, Ferrenberg and Landau (CFL) 0, who 
studied the 8-state random-bond Potts model with a bi- 
modal self-dual distribution. On the other hand the ex- 
perimental evidence has been found that, for the order- 
disorder phase transitions of absorbed atomic layers on 
two-dimensions, the critical exponents are changed from 
the original universality class of the 4-state Potts model 
on the addition of disorders [^|, gj. However, no modifica- 
tion is found when the pure system belongs to the Ising 
universality class fTpf . Recently such disordered systems 
are extensively studied by examining how a phase tran- 
sition is modified by quenched disorder coupling to the 
local energy density , where use of intensive MC sim- 
ulations is often helpful @ g ||, || ||. 

The two-dimensional (2D) q-state random-bond Potts 
ferromagnet (RBPF) is an interesting framework in the 
MC research to study the influence of impurity on pure 



systems. For q > 2 such randomness acts as a relevant 
perturbation, and for q > 4 it even changes the nature 
of the transition from first to second order. In Table 1 
we list the magnetic scaling index (3/v of the 2D 8-state 
RBPF obtained by different groups. Here a disorder am- 
plitude r, defined by the ratio of the strong to weak cou- 
pling (distributed according to the bimodal distribution) 
in the range r =2-20 appears to be adapted to a numer- 
ical analysis and gives a good estimate of the disordered 
fixed point exponents [§, [17|] . A recent work of Olson and 
Young |l6j used a specially continuous self-dual probabil- 
ity distribution of the disoreded bonds 
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, and performed a MC 



study of multiscaling properties of the correlation func- 
tions for several values of q. Their results are very in- 
teresting to examine the universality class of the RBPF. 
Cardy and Jacobsen studied the RBPF based on the 
connectivity transfer matrix (TM) formalism [^9|, and 
their estimates of the critical exponents lead to a contin- 
uous variation of (3/v with q, which is in sharp disagree- 
ment with the MC results by CFL(see Table 1). 

In this paper, we present a MC study by short-time 
dynamics (STD) to verify the dynamic scaling features 
of the RBPF and estimate the critical exponents for a 
trinary random-bond Potts model on a 2D triangular 
lattice. It is well known that this model on the 2D 
triangular lattice belongs to the same universality as 
the one on the 2D square lattice in equilibrium. How- 
ever, a difference exists that the transition temperature 
is changed because the systems on triangular lattices 
lost the original (square-square) self-duality relation. In- 
stead, a honeycomb-triangular duality (or star-triangular 
duality) is satisfied p2| . Therefore the multi-disorder 
amplitudes r.; can be introduced to study the dynamic 
universalit y o f the RBPF on such lattices by the STD 
approach (231 E4j. In particular we will investigate the 



2 



Authors 



P/v Technique 



CFL 

Chatelain and Berchcfll5 ] 



2,10 0.118(2) MC 
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0.1505(3) 


TM 
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TABLE I: Magnetic exponent fj/v estimated by different 
groups for the 2D 8-state RBPF. 



dynamic behavior of critical scaling affected by introduc- 
ing quenched randomness, and study the dependence of 
critical exponents on the randomness strength in order 
to clarify the crossover behavior (pure < — ► random fixed 
point < — ► percolation-like limit) ]l7[. 



II. THE MODEL 

The Hamiltonian of the g-state Potts model on a 2D 
triangular lattice with quenched random interactions can 
be written as 
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where the spin <7j, defined on lattice site i, takes the 
values 1, • • • q, /3 = l/ksT is the inverse temperature, 5 
the Kronecker delta function and the sum is over all the 
nearest-neighbor (NN) interactions (bonds) < i,j > on 
the lattice with a size N — L 2 . There are 3iV NN bonds 
(six for each and shared by two spins) and they can be 



J 



grouped into three classes {Lq, Li, L 2 } corresponding to 
the interactions on them. The dimensionless couplings 
Kij can be selected from three positive (ferromagnetic) 
values Kq, K\ = r±Ko and K 2 = r 2 Ko according to a 
probability distribution, 



P(K) = (1 - Pl - p 2 )8(K - K Q ) + Pl S(K -Kx)+ P2 S(K - K 2 ). 

I 



(3) 



When pi — p 2 — 1/3 it is a trinary random-bond system. 
By denoting F = (e K - - 1), Fx = (e< - 1) and F 3 = 
(e K <= — 1), where K c , K' c = rxK c and K" = r 2 K c are the 
critical values of Ko, Kx and K 2 respectively, the critical 
point can be determined by a self-dual relation p3, E5| , 



FqFxF 2 + -Fbfr + FqF 2 + FxF 2 — q 



(4) 



for the trinary system with the given state parameter q. 

The strong to weak coupling ratio = Ki/Ko is called 
disorder amplitude. The values r\ = or r 2 = corre- 
spond to the diluted case ^6|, and rx = r 2 = 1 to the 
pure case where the phase transitions are first-order for 
q > 4. With the presence of quenched random-bonds, 
however, the second-order phase transitions are induced 
for any q-state of the Potts model and the critical points 
can be calculated according to Eq.(|J). In our simulations 
we chose q = 8 which is known to have a strong first- 
order phase transition in the pure cases, in hope that 
we would find second-order phase transitions induced by 
the quenched random disorders to demonstrate the in- 
fluence of quenched impurities on the first- order phase 



transitions. The strength of disorder is realized by the 
disorder amplitudes {ri,r 2 } as chosen in Table 2. There 
are three types of such disordered systems for the tri- 
nary distribution: (I) one-third bonds chosen randomly 
are strongly coupled with {ri=l, r 2 > 1}, (II) two- 
third bonds strongly coupled with different amplitudes 
T% > rx > 1 and (III) two-third bonds with the same am- 
plitude r\=r 2 > 1. Actually the systems III, 2/3 bonds 
being strongly coupled with {rx=r 2 > 1}, are related 
to those I, 1/3 bonds being strongly coupled, because 
the former can be transferred to the latter by taking 1/3 
bonds to be weakly coupled with {r*i=l, r 2 = l/r 2 < 1}. 
The study will be concentrated on the important ques- 
tion whether there exists an Ising-like universality class 
for systems with multi-quenched randomness |l3l p20[ . 



III. SHORT-TIME DYNAMICS 

For a long time, it was believed that universality and 
scaling relations can be found only in equilibrium or in 
the long-time regime. In Ref . [p3[ , however it was dis- 
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{ri,r 2 } 


T c 


{ri,r 2 } 


T c 


{1,5} 


.39903934... 


{2, 5} 


.33624310... 


{5, 5} 


.19047260... 


{1,8} 


.29663628... 


{5, 8} 


.19047260... 


{8, 8} 


.15574523... 


{1,10} 


25548102... 


{8,10} 


13987343... 


{10,10} 


.12638846... 


{1,12} 


22533080... 


{2,12} 


19784059... 


{12,12} 


.11306096... 



TABLE II: The critical temperature as a function of disorder 
amplitudes {ri, r 2 } for the trinary RBMF. (I) 1/3 bonds are 
strongly coupled with r 2 > 1, (II) 2/3 bonds are strongly 
coupled with different amplitudes r 2 > n > 1 and (III) with 
the same amplitude r\=r 2 > 1. 



covered that for an O(N) vector magnetic system in the 
states with a very high temperature T » T Cl when it 
is suddenly quenched to the critical temperature T c and 
evolves according to a dynamics of model A, a univer- 
sal dynamic scaling behavior emerges already within the 
short-time regime, 

M<*'(t, r, L, too) = 6" fc ^"M ^ fc 1 ^ b' 1 L,b x °m ), 

(5) 

where M^ k ' is the fcth moment of magnetization, t,r — 
(T — T c )/T c , L and b are time, reduced temperature, lin- 
ear size of the lattice and scaling factor respectively. (3 
and v are the well known static critical exponents. The 
quantity xq, a new independent exponent, is the scaling 
dimension of the initial magnetization Too- This dynamic 
scaling form is generalized from finite size scaling in equi- 
librium. 

Up to present, the short-time dynamic MC simulations 
have been successfully performed to estimate the critical 
temperatures T c and the critical exponents 9, [3, v and 
the dynamic exponent z for the various models with sec- 
ond order phase transitions |24|] . This approach has also 
been extensively applied for systems with disorders, as 
the FFXY and spin glass, to estimate both dynamic and 
static critical exponents 



28 , and for systems on the 



2D triangular lattices to study the dynamic scaling uni- 
versality p9| |. The method is also efficient to study the 
deterministic dynamics |]§o| . 

We carry out simulations in an analogous strategy de- 
scribed in Ref. |^l| where the magnetic observables to be 
measured for the RBPF have been well defined. In this 
work, however, an additional step is introduced to verify 
that the second-order phase transitions really are induced 
when quenched randomness is added. Then we will be 
particularly concerned to demonstrate the evidence for 
the dynamic universality of the trinary RBPF. Therefore 
our MC simulations are mainly performed for the time 
evolution of the magnetization M(t) from the disordered 
initial states with small magnetizations too > 0, and the 
completely ordered state with too = 1. In the former, 
the magnetization M(t) undergoes an initial increase at 
the critical point K c for t > t m iej 



For the latter, with Too 
decay behavior Q 

M(t) ~ r c \ 



1, M(t) shows a power-law 



ci = /3/vz, 



(7) 



at the critical point. Furthermore the Binder cumulant 
U(t,L) = M^(t,L)/M 2 (t,L) shows, for m = 1 at crit- 
icality, a similar power-law behavior on a large enough 
lattice, 



U(t,L)~t Cu , c u = d/z. 
Here the magnetization M(t) is defined by 

.qM (t)-N 
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M(t) ~ m t 8 
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AI = max(Mi, M2, • • • , M q ) with Mi being the number 
of spins in the ith state. = L 2 is the total number of 
spins on the lattice. < • • • > denotes thermal averages 
over the initial states and/or random number sequences, 
and [• • •] the disorder averages over quenched randomness 
distributions. The time unit t is defined as a MC sweep 
over all spins on the lattice. 

The timc-dcpcndcnt susceptibility (the second moment 
of the magnetization) is defined as usual 

M< 2 >(i) = i [(< M 2 {t) > - < M{t) > 2 )] . (10) 

It has played an important role in equilibrium to deter- 
mine the critical exponents j/v and (3/v @]. 



IV. RESULTS 

In our simulations, up to 60,000 MC processes are 
taken for MC averages (300 samples chosen as disorder 
averages in a given distribution and 100-200 initial con- 
figurations as thermal averages for each disorder config- 
uration realization). Statistical errors are simply esti- 
mated by performing three groups of averages with dif- 
ferent random number sequences as well as independent 
initial configurations. To minimize the number of bond 
configurations needed for the disorder averaging, we con- 
fine our study to the bond distributions in which there 
are the same number of disordered (stronger) NN bonds 
for each of the three groups {Lo, L±, L2}. This procedure 
is reliable and should reduce the variation between dif- 
ferent disorder configurations with no loss of generality. 
Periodic boundary conditions are applied to the N = L 2 
lattice. We use L up to 128 and adopt the heat-bath 
algorithm. 

The data for the evolution of M(t) from too > and 
too = 1 exhibit that the power-law behavior is satisfied 
at the same point K c for both the disordered and ordered 
initial states. This fact gives an evidence that a second- 
order phase transition is induced at K c instead of a first- 
order one: In Ref.[|4| it was argued that there are two 
quasi-critical points if a first-order transition happens, 
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FIG. 1: M(t) versus MC time plotted on a log-log scale for 
{n=l,r2=8} from the initial states mo ~ 0.01, performed at 
and near the critical point K c and = (1 ± 0.012)if c . 



FIG. 3: M(t) versus MC time plotted on a log-log scale for 
{ri=5,r2=8} from random initial states with magnetizations 
m =0. 04,0. 02,0. 01 and with m ~0ona 128 2 lattice. 
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FIG. 2: The same as Fig.l but from the ordered states mo = FIG. 4: M(t) versus MC time plotted on a log-log scale for 
1. {n = l,f"2=8}, {ri=l,r2 = 10} and {ri = l,r2=12} from the ini- 

tial states mo = 1 on a 128 2 lattice. 



■^■random f or evo i u tions from random initial states and 
K order from or d er ed initial states, and K r c andom < K c < 
K° rder is satisfied. We plot in Fig.l and Fig. 2 the curves 
M(t) for disorder amplitudes {n=l, r 2 = 8} on a 64 2 
lattice. They exhibit obviously a power-law relaxation at 
the same critical point for both the mo ^ and m = 1 
states. On the other hand the curves deviate from the 
power-law behavior for a small but nonzero r described 
by M(i,r) = m t 6 'M (t 1 ' ll ' z V) or M(t, r) = t- Ci F(t x l vz T) 
with respect to the disordered or ordered initial states 
respectively. This modification has been initially used to 
determine the critical temperatures ]32] , |33| ]. The same 
features, as plotted in Figs. 1 and 2, exist for all other 
realizations of the randomness listed in Table 2. 



Then we perform simulations systematically at the ex- 
act critical points K c {ri) as a function of {^1^2} for 
magnetization evolutions starting with mo ^ and 
mo = 1. The disordered initial states are prepared from 
random states with small magnetizations mo = 0.04 - 
0.01. The curves M(t) for the disordered configurations 
with r\ = 5, T2 = 8 on a 128 2 lattice are displayed in 
Fig. 3 with a double-logarithmic scale. They exhibit the 
power-law behavior of Eq.(^) for t > t m i C , which depends 
on mo. Thus the exponent 9 can be estimated with the 
standard least square fitting algorithm, from the slopes 
of the curves for t > t m i C , and by an extrapolation to 
the mo — > limit. The values are included in Table 3. 
Next we start from mo = 1 to observe the power-law de- 
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{Vi,r a } 
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Cl c u 


z &/v 


I: {1,5} 
{ 1,8} 

{ 1,10} 

{ 1,12} 


.238(3) 
.206(2) 
.187(2) 
.180(2) 


.0780(8) .789(6) 
.0654(7) .611(6) 
.0597(7) .599(6) 
.0559(7) .554(6) 


2.53(4) .202(5) 
3.28(5) .212(5) 
3.34(5) .200(5) 
3.61(5) .202(5) 


II: {2,5} 
{5, 8} 
{ 8,10} 


.268(3) 
.283(3) 
.290(3) 


.0864(9) .889(6) 
.0849(9) .893(7) 
.0830(9) .897(7) 


2.25(4) .194(5) 
2.24(4) .190(4) 
2.23(4) .185(4) 


III: { 5, 5} 
{8, 8} 
{10,10} 


.305(4) 
.296(3) 
.293(3) 


.0784(8) .959(7) 
.0811(9) .907(8) 
.0815(9) .898(8) 


2.09(3) .164(4) 
2.21(4) .179(4) 
2.23(4) .182(4) 



TABLE III: The values of indices (8,ci,c u ), and results of 
exponents (z, /3/v) with respect to the disorder amplitudes 

{n,r 2 }. 



cay of Eq.(|7|) at the critical points K c (ri). Some curves 
M{t) are plotted in Fig. 4, which appear to have also a 
nice power-law behavior, and the slopes are then used to 
estimate the index c\ = (3/vz (for results see Table 3). 

Now we turn to the dynamic scaling behavior of the 
time-dependent Binder cumulant from mo = 1 to esti- 
mate the scaling index c„ = d/z. Unlike the relaxation 
from the disordered state, the fluctuations involved in the 
ordered initial states (e.g., Ui = 1, i = l,N) are much 
smaller. Therefore, the measurements of the critical ex- 
ponents based on Eq. (|7j) and Eq.(||) are better in quality 
than those from disordered states on Eq.@. So we take 
about 100 MC sweeps for random number sequences to 
be the thermal averages for the completely ordered ini- 
tial state. In Figs. 5 the curves U(t) on a 128 2 lattice 
are presented and they exhibit a power-law behavior at 
the critical points as a function of disorder amplitudes 
{n, r 2 }. Finally, based on the results of c u = d/z and 
ci = (3/vz, the critical exponents z and (3/v can be cal- 
culated consistently and their results are summarized in 
Table 3. 



V. SUMMARY AND CONCLUSION 

We have, for the first time, introduced the trinary dis- 
tribution for the 2D RBPF with multi-disorder ampli- 
tudes to study its critical behavior by the STD method. 
In the MC simulation the dynamic scaling behavior is 
verified and the power-law behavior is used to estimate 
both the dynamic and static exponents as a function of 
the disorder amplitudes {r\, r 2 }. It is found that the val- 
ues of exponents 6 and (3/v vary continuously with the 
disorder amplitudes, and they violate the Ising-like uni- 
versality. Furthermore the values of dynamic exponent 
z are larger than those for systems without disorders 
and they increase with the strength of disorder ampli- 
tudes {ri,f2} for the systems I: r± = 1, r% > 1 and III: 



ri=?*2 > 1 in Table 3. 




MC time 



FIG. 5: U(t) versus MC time on a log-log scale as a function of 
the disorder amplitudes {ri, T2\ from the initial state mo=l, 
used to estimate the index c u . 
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FIG. 6: The U(t) curves for the disorder amplitudes {n=5, 
V2=8}, {ri=8, r2=10} and {ri=10, r2=10} are almost paral- 
lel each other, but it dose not for the {ri=l, r2=5} curve. 

On the other hand, for the systems II where 2/3 bonds 
are strongly coupled with different amplitudes 1 < n < 
r2, they seem to have been located at a "random" regime 
where their values of the critical exponents (3/v and z are 
nearly independent of disorder amplitudes within the er- 
ror bars, z ~ 2.23 and (3/v ~ 1.85. By comparing these 
results to z=2.23(4) and /3/i/=0.182(5) for the ri=r 2 = 10 
in III, we could argue that, after r\ = r 2 > 10 (equiva- 
lent to ri=l, r 2 < 0.1 in I) it will pass the "crossover" 
region to the random regime from other type, I or III, 
of the disordered systems, so that z and (3/v will be 
nearly constants. In Fig. 6 the U(t) curves for the dis- 
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order amplitudes {ri — 5,r2 = 8}, {ri = 8,r2 = 10} and 
{rj = 10, r 2 — 10} are almost parallel each other, which 
shows that the values of exponent (3/v are nearly same 
for these disordered systems. 

In conclusion, our simulation verifies that second or- 
der phase transitions are induced. Then the results show 
evidence that the dynamic universality class of the tri- 
nary RBPF should not belong to that of the Ising model, 
as inferred by Jacobsen and Cardy |19| ]. From the work 
we find it rather encouraging to apply the short-time dy- 
namic MC to simulate the scaling and critical dynamics 



of disordered spin systems. We will pursue this field to 
explore the logarithmical slow dynamics |^|, |35| . 

Acknowledgement: We are grateful to the helpful dis- 
cussions with Q. Wang. H.P.Y would like to thank the 
Heinrich-Hertz-Stiftung for fellowship and acknowledge 
the hospitality of Universitat Siegen. The work sup- 
ported in part by the Deutsche Forschungsgemeinschaft, 
DFG Schu 95/9-3 and by the NNSF of China, 19975041 
and 10074055. 



J. Cardy, in Scaling and Renormalization in Statisti- 
cal Physics, (Cambridge University Press, Cambridge, [19 
1996). 

J. Cardy, Physica A263, 215 (1999), and references there [20 
in. 

A.B. Harris, J. Phys. C7, 1671 (1974). [21 
Y. Imry and M. Wortis, Phys. Rev. B19, 3580 (1979). [22 
K, Hui and A.N. Berker, Phys. Rev. Lett. 62, 2507 
(1989). [23 
M. Aizenman and J. Wehr, Phys. Rev. Lett. 62, 2503 
(1989). [24 
S. Chen, A.M. Ferrenberg and D.P. Landau, Phys. Rev. [25 
Lett. 69, 1213 (1992); Phys. Rev. E52, 1377 (1995). 
L. Schwenger, K. Budde, C. Voges and H. Pfhiir, Phys. [26 
Rev. Lett. 73, 296 (1994). 
C. Voges and H. Pfniir, Phys. Rev. B57, 3345 (1998). [27 
Ch.V. Mohan, H. Kronmiiller and M. Kelsch, Phys. Rev. 
B57, 2701 (1998). [28 
R. Folk, Yu. Holovatch and T. Yavors'kii, Phys. Rev. 
B61, 15114 (2000). [29 
S. Wiseman and E. Domany, Phys. Rev. E51, 3074 
(1995). [30 
M. Kardar, A.L. Stella, G. Sartoni and B. Derrida, Phys. 
Rev. E52, R1269 (1995). [31 
F. Yasar, Y. Giinduc and T. Celik, Phys. Rev. E58, 4210 [32 
(1998). [33 
C. Chatelain and B. Berche, Phys. Rev. E58, R6899 
(1998). [34 
T. Olson and A. P. Young, Phys. Rey . B60, 3428 (1999). [35 
M. Picco, e-print |zond-mat/9802092j . 
J. Cardy and J.L. Jacobsen, Phys. Rev. lett. 79, 4063 



(1997) . 

J.L. Jacobsen and J. Cardy, Nucl Phys. B515[FS], 701 

(1998) . 

C. Chatelain and B. Berche, Phys. Rev. E60, 3853 

(1999) . 

H.-P. Ying and K. Harada, Phys. Rev. E62, 174 (2000). 
R.J. Baxter, in Exactly solved models in Statistical Me- 
chanics (Academic Press, London, 1982). 
H.K. Janssen, B. Schaub and B. Schmittmann, Z. Phys. 
B73, 539 (1989). 

B. Zheng, Int. J. Mod. Phys. B12, 1419 (1998). 

D. Kim and R.I. Joseph, J. Phys. C: Solid State Phys. 7, 
L167 (1974). 

M.I. Marque and J.A. Gonzalo, Phys. Rev. E60, 2394 
(1999); E62, 191 (2000). 

H.J. Luo, L. Schiilke and B. Zheng, Phys. Rev. Lett. 81, 
180 (1998). 

H.J. Luo, L. Schiilke and B. Zheng, Mod. Phys. Lett. 
B13, 417 (1999). 

H.-P. Ying, J. Wang, J.B. Zhang, M. Jiang and J. Hu, 
Physica A, to be published. 

B. Zheng, M. Schulz and S. Trimper, Phys. Rev. Lett. 
82, 1891 (1999). 

Nobuyasu Ito, Physica A 192, 604 (1993). 
L. Schiilke and B. Zheng, Phys. Lett. A204, 295 (1995). 
Jaster, J. Mainville, L. Schiilke and B. Zheng, J. Phys. 
A: Math. Gen. 32, 1395 (1999). 

L. Schiilke and B. Zheng, Phys. Rev. E62, 7482 (2000). 
H. P. Ying, B. Zheng, Y. Yu and S. Trimper, Phys. Rev. 
E63, R035101 (2001). 



